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A rigorous formulation for the scattering of surface plasmon polaritons (SPP) from a one-dimensional 
surface defect of any shape that yields the electromagnetic field in the vacuum half-space above the 
vacuum-metal interface is developed by the use of an impedance boundary condition. The electric 
and magnetic near fields, the angular distribution of the far-field radiation into vacuum due to SPP- 
photon coupling, and the SPP reflection and transmission coefficients are calculated by numerically 
solving the fc-space integral equation upon which the formulation is based. In particular, we consider 
Gaussian-shaped defects (either protuberances or indentations) and study the dependence of the 
above mentioned physical quantities on their 1/e half- width a and height h. SPP reflection is 
signiflcant for narrow defects (a ~ A/5, for either protuberances or indentations, where A is the 
wavelength of the SPP); maximum reflection {plasmon mirrors) is achieved for a ~ A/10. For 
increasing defect widths, protuberances and indentations behave differently. The former give rise to 
a monotonic increase of radiation at the expense of SPP transmission for increasing defect half-width. 
However, indentations exhibit a significant increase of radiation (decrease of SPP transmission) for 
half-widths of the order of or smaller than the wavelength, but tend to total SPP transmission in 
an oscillatory manner upon further increasing the half-width. Both the position of the maximum 
radiation and the oscillation period depend on the defect height, which in all other cases only affects 
the process quantitatively. Light- emitters might thus be associated with either wide indentations, 
or protuberances with widths that are of the order of or smaller than the wavelength. 

PACS numbers; 73.20.Mf, 61.16.Ch, 78.66.Bz, 42. 25. Fx 



I. INTRODUCTION 

In this paper, we study the scattering of surface plas- 
mon polaritons (SPP) by surface defects. SPP are 
p-polarized electromagnetic (EM) waves bound to a 
dielectric-metal interface and caused by thcr^urface oscil- 
lations of the electron plasma of the metahlil They prop- 
agate along the metal interface a distance of the order 
of the SPP mean free path (ranging from microns in 
the visible to millimeters in the infrared, of course de- 
pending also on the metal being considered), undergoing 
scattering processes due to surface roughness. This con- 
stitutes a classical problem of fundamental interest not 
only in the case of individual defects (cf. Ref. || and ref- 
erences therein), but also forppfiLniodically or randomly 
(or both) distributed defects. Qm'El Furthermore, it is ob- 
viously crucial in any light scattering problem involv- 
ing rough metal surfaces where roughness-induced exci- 
tation of SPP occurs. This has been sjxjJicitly shown 
in connection ||W.itb. either single defectsBuQ or random 
corrugation,E3Eilll3'E3 the latter configuration being rel- 
evant to the phenomenon of (SPP mediated) enhanced 
backscattering of light. In addition to that, light-SPP 
coupling plays a central role in other phenomena such as 



anomalaus, transmission through metal slabs withjitole 
arrays sucf ace-enhanced Raman scattering 
or biosensingJl^ 

In recerd— years, the advent of near-field optical 
microscopyEj has opened up the possibility to study ex- 
perimentally SPP in a direct manner. Among the var- 
ious configuratians developed, photon scanning tunnel- 
ing microscopycJ (PSTM), basically exploiting SPP ex- 
citation in the attenuated total reflection arranger nent . 
has made it possible to probe the SPP struct rmjHS'tJ 
localized SPP on randomly rough surfaces JHj and 
SPP resonaiu:n£i^ in fractal colloid clustersEj and sin- 
gle particles .oo Moreover, PSTM images have been 
obtained by surface-enhanced Raman scattering prob^ 
ing single molecules adsorbed on single nanoparticles.c3 
PSTM in combination with direct-write lithography has 
made it jppssible to create sub-micron defects on metal 
surfaces. E3 

Particularly relevant to the present work are the re- 
cent experiipept^ studies on SPP scattering by sur- 
face defects. E3E3'El These studies have shown evidence 
of drastically distinct scattering properties depending 
on the defect size. Specifically, surface defects favoring 
SPP reflection and light coupling, called SPP mirrors 



and flashlights Jsi respectively, have been ^described, as 
well as SPP microlenses and microcavitiesj^ SPP Bloch 
waves have also been imaged in periodic arrays of surface 
defectsEil Interestingly, the possibilities of artificially cre- 
ating micro-optical components for SPP have also been 
noted in these studies. Much has to be done, however, 
from the theoretical standpoint. Quite recently, calcu- 
lations for circularly symmetric defects have successfully 
accounted for tha-*)eculiar azimuthal dependence of the 
radiated pattern;E3 in addition, such calculations have 
been used to retrieve the surface profileEj In the case 
of one-dimensional surface defects, preliminary calcula- 
tions have focused on the optimization of the defecLsize 
to obtain SPP mirrors and so-called light-emitters.o In 
this regard, it is our purpose to address in detail the SPP 
scattering by one-dimensional surface defects, including 
near-field and far-field calculations (along with energy 
balance) and their dependence on defect size parameters. 
Thus we expect not only to shed light on the experimen- 
tal works mentioned above, but also to find and predict 
related effects. 

The physical system we consider here is a planar one- 
dimensional metal surface with a one-dimensional de- 
fect. The surface corrugation is modeled by using a local 
impedance boundary condition (IBC) on a fiat surface. 
The connection between surface impedance and reial sur- 
face corrugation has been recently demonstrated,EJ and 
its validity to give accurate quantitative results has been 
shown in numerical calculations of grating-induced SPP- 
photon coupling.Ll A scattering-theoretic formulation of 
the interaction of an SPP with the surface roughness 
is developed by imposing the IBC on the amplitude of 
the magnetic field in the vacuum region in the form of a 
Rayleigh expansion. Upon solving the resulting integral 
equation for the scattering amplitude, the magnetic field 
at any point in the vacuum half-space can be calculated. 
We will focus on the far field angular distribution and 
the surface field amplitudes to determine, respectively, 
the total radiated energy S, and the SPP refiection Rsp 
and transmission Tsp coefficients. By numerical simula- 
tion calculations, these quantities are computed. 

The paper is organized as follows. The theoretical for- 
mulation is derived in Sec. and some details pertaining 
to the numerical procedure are given in the Appendix. In 
Sec. [II, we show the results obtained for a single Gaus- 



sian defect and the influence of defect width and height. 
Finally, Section |l^ summarizes the conclusions drawn 
from this research. 



II. THEORY 



A. Scattering Equations 



We study the scattering of a p-polarized SPP of fre- 
quency Lu propagating along a flat vacuum-metal inter- 
face (x3 = 0) by a one-dimensional obstacle (constant 



along the a;2-axis, see Fig. |^). Under these circumstances, 
the three-dimensional electromagnetic problem can be 
cast into a two-dimensional scalar problem in such a way 
that the single, nonzero component of the magnetic field 
amplitude H2{xi, X3) is the solution of the corresponding 
two-dimensional Helmholtz equation in the upper (vac- 
uum) and lower (metal) half spaces. The magnetic field 
in vacuum is assumed to be the sum of an incoming SPP 
and a scattered field as follows 

H^{xi,X3) = ex.p[ik{uj)xi - /3o(a;)x3] 



dq 
2^ 



R{q,LLi)eyi-p[iqxi+iao{q,uj)x3], (1) 
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Note that the expressions for the SPP wavevector com- 
ponents k{uj) and Po{uj) in vacuum apply in the limit 
e(a;)| 3> 1. This stems from the fact that the continuity 
conditions across the interface are mapped onto a local 
IBC on the planar surface X3 = in the form 
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0x3 
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where the superscript > indicates the vacuum region, 
— (w/c)[— e(cij)]~-^/^s(xi) is the contribution to the sur- 
face impedance associated with the obstacle, and e(a;) is 
the isotropic, frequency-dependent dielectric function of 
the metal. The IBC has been widely used in the past 
to model the vacuum-metal interface qualitatively, es- 
pecially in the infrared region of the optical spectrum. 
Furthermore, it has been recently proven to be quantita- 
tively accurate in calculations of grating-induced photon- 
SPP couplingj by using the connection between surface 
impedance and real corrugation demonstrated in Ref. |3^. 

In order to calculate the scattering amplitude R{q,uj), 
we substitute Eq. ([l]) into Eq. (^), and obtain the inte- 
gral equation R{q, w) satisfies, 

Riq,uj) = Go{q,Lu)V(q\k{Lo)) 

+G^{q,u:) r ^V{q\p)R{p,u), (6) 
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where 

Go{q,uj) 



5Hao(Q,a;)+z(u;/c)[-eM]i/2 



(7) 



is the Green's function of the SPP on the unperturbed 
surface [s{xi) — 0]. We have also introduced the scatter- 
ing potential 



V{q\p)=(3o{^)s(q~p), 



with 



m = / 



dxi e 



(8) 



(9) 



to simplify the notation. Equation (^) can be rewritten 
in a more convenient manner by substituting 



R{q,uj) = Go{q,Lj)T{q,uj), 
into it, so that 

T{q,u;)^V(q\k{L,)) 



(10) 



dp 
,2^ 



V{q\p)Go{p,Lo)T{p,u;). (11) 



Equation (^, along with Eqs. (|^) and (pX|), is the basis 
of our theoretical formulation. 

In solving Eq. (p]), it is very important how we deal 
with the poles appearing in the Green's function (|^). 
First, we rewrite the latter in the form 



Go(q,w) = C{q,ui) 
with 

C{q,u; 



1 



1 



q — k[uj) q + k{uj) 

^ e{Lu)ao{q,u;) - z(tj/c)[-e(^)]^/^ 
2ie(Lu)k{uj) 



(12) 



(13) 



We now assume that the metal dielectric function is given 
by Drude's expression 



(14) 



where ujp is the plasma frequency, in the absence of ab- 
sorption losses. Therefore, in light of Eq. (^, we have to 
take the limit fc^(aj) — > in Eq. (|^), to obtain 



Go((?,w)= C{q,uj) 



q~k^{uj) 



k^iuj) 



+ Tri[S{q - fc^(w)) + Siq + k^'iuj))] 



(15) 



The first two terms on the right-hand side of Eq. ( p^ ) 
have meanings in the Cauchy's principal value sense, 
whereas the last two terms are delta functions. Once 
Eq. ( |Tl| ) is solved for T{q,oj) [we will see below how to 
do so numerically with the help of Eq. (p^)], we pro- 
ceed to calculate the electric and magnetic near fields, 
the SPP-photon coupling, and the SPP reflection and 
transmission coefficients in the following manner. 



B. Near Field 

The magnetic field at any point in the vacuum half- 
space can be straightforwardly calculated from Eq. (|^), 
upon recalhng Eq. (|lo|), which relates T{q,uj) with the 
scattering amplitude R{q, lo). Then the electric field com- 
ponents in vacuum are easily written also as functions of 
i?(q, w) by means of a Maxwell curl equation as follows: 



-Ef (xi,X3) = -i(3o{uj) ey.^[ik{uj)xi -/3o(w)x3] 

LO 



dq 



aQ{q,uj)R{q,uj) 



^ J -co 27r 
X e:}q3[iqxi + iao{q,uj)x3] (16a) 
E>{xuX3) = (16b) 

E^{xi,X3) = k{uj) ex.p[ik{uj)xi - l3o{uj)x3] 



LO 

C dq 



LO 2tt 



qR{q,io) 



X eyipliqxi + iao{q,uj)x3]. (16c) 

The time-averaged Poynting vector thus reads: 

< S >= (E X H*) = {-E3H;,0, EiH;) , (17) 

where 5ft denotes the real part and the asterisk the com- 
plex conjugate. 

C. Radiated energy 

The total power carried away from the surface in the 
form of volume electromagnetic waves propagating in the 
vacuum region above it, per unit length of the system 
along the X2-axis is 



P,r. - / dxi < Slf"^ > 



oo 

2 fi^/c 



Sttuj 



I p-ao{q,Lo) |i?((j,w)|^ 



(18) 



Note that only the scattered field contribution to the X3- 
component of the time-averaged Poynting vector is used. 
Equation ( p^ ) must be normalized by the power carried 
by the incident SPP per unit length along the a;2-axis 



P,...r = I dx3 < S['"^'> >= 



where < S'^"^-' > is the xi-component of the time- 
averaged Poynting vector of the incident SPP. Then the 
total, normalized scattered power S is given by 



c^kjio) 
167r[j/3o(w) ' 



(19) 



(20) 
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where 



dR 1 I3q{uj) 



2tt 2k{uo) 



R[q= —su 
\ c 



(21) 



is the differential reflection coefficient (DRC), namely, 
the fraction of the energy of the incident SPP that is scat- 
tered into an angular region of width dOs about the scat- 
tering direction 6s where the scattering angle 9s is mea- 
sured clockwise with respect to the xa-axis (see Fig. ||). 



D. Reflection and transmission coefficients 



points on the mesh, as required by Eq. ([231). In ad- 
dition, the discretization is not regular: the density of 
q-points around the poles aX q — ±k^{u!) is considerably 
larger (Aq w 10~^uj/c) than it is either in the radiative 
region \q\ < oj/c or in the non-radiative region away from 
the poles {Aq « 10~^w/c). The number TV of g-points 
needed in the numerical procedure depends not only on 
the accuracy required to sample the pole regions, but 
also on the explicit form of the obstacle, which enters in 
the calculation through its Fourier-transform in Eq. 
Throughout this work, typically N = 2600, except for 
the larger defects for which up to = 4000 points are 
employed. The convergence of the numerical results with 
increasing N has been checked in the most unfavorable 
cases. 



In order to evaluate the amplitude of the reflected and 
transmitted SPP, we study the behavior of 7J^(xi,a;3), 
Eq. (0), with the help of Eqs. (10) and (15), on the 
surface — 0. At this point, great care has to be taken 
when calculating the contribution to the scattered field in 
Eq. (2) from the Cauchy principal value integrals arising 
from the first two terms on the right-hand side of Eq. 
(15). We assume that the obstacle has a finite extent 
and is centered about Xi = 0. If we focus on the regions 
Xi -C and xi ^ far from the obstacle, it can be shown, 
by working out the contributions from those integrals ia, 
the complex q-space with the help of Cauchy's theorem,E3 
that the magnetic field is given by: 

iJ^(xi,X3 = 0) = exp[ik^'{u!)xi] 

+r(tj) exp[-ifc-"(w)xi],xi < (22a) 
= t{Lu) exp[ifc^(w)xi], xi > 0, (22b) 

where the amplitudes of the reflected and transmitted 
SPP, r{uj) and t{uj), respectively, are 



r(cj) = iT{ - fc^(w),w) C{ ~ k^{w),uo) 
t{uj) = 1 + iT{k'^{uj),oj) C{k'^{u),oj). 



(23a) 
(23b) 



Equations (|2|) manifest the fact that, away from the ob- 
stacle, only the incident and reflected SPP (on the left- 
hand side, see Fig. 1) and the transmitted SPP (on the 
right-hand side) propagate along the interface. The cor- 
responding reflection and transmission coefficients are 



^H-k(^)P 

T{uj) = \ t{uj)\^ . 



(24a) 
(24b) 



E. Numerical Calculations 



The integral equation ( p_lD is numerically solved by 
converting it into a matrix equation through a quadra- 
ture scheme. The details are given in the Appendix. 
It should be pointed out that the discretization g-mesh 
is chosen in such a way that q = ±k^{u!) are always 



III. RESULTS AND DISCUSSION 

Note that up to now no restrictions have been imposed 
on the shape of the obstacle apart from its having a finite 
extent along the xi-axis. Its surface impedance function 
s(xi) is connected to-jthe actual surface profile defined by 
X3 = /(xi) throughEa 



s(x) 



1 - e(^) 
d{u!)e{uj) 



[l^d^iu)D']'/^fix^) + Oif), (25) 



where d{iL!) — (c/ e{Lij)]^^ is the optical skin depth, 
and D = d/dxi. In the case of small skin depths and 
surface slopes (dD)^ <C 1, the square root term on the 
rhs of Eq. (|2^) can be expanded as 



211/2 



[1 - [dDY] 
1.1.3. 



1 



1 



2 

(2n - 3) 



2.4.6. . .2r7, 



{dOr^... 

{dDf' + 0{{dDf''+^). (26) 



Then the Fourier transform of the surface impedance 
function, which is needed in the calculation [cf. Eq. (||)], 
is related to the Fourier transform f{Q) of the surface 
profile function through 



m = - 



1 - 

d{uj)e(uj) 



l-\[^id{u)QY 



- -[-^d{u)Q]^ + 0{[-id{u)Qf)J f{Q). (27) 

In what follows, we will restrict the analysis to a Gaus- 
sian defect of 1/e half- width a and height h: 



/(xi) = ft-exp [-x\/a^) 



(28) 



In addition, unless otherwise stated, we retain in Eq. 
(p6|), and thus in Eq. (|27|), only the zeroth order term 
in the expansion in powers of [—id{ijj)Q]'^ , as implicitly 
done in Ref. Therefore the function s(Q) we will use 
in our calculations is 
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s(Q) =^i/2soaexph(aQ)V4], 



with 



So = 



1 - e{uj) h 
t{uj) d{uj) 



(29a) 



(29b) 



transmission coefficients are thus calculated from Eqs. 
(H) and (I2I 



It should be emphasized that the approximation involved 
in retaining only the lowest order term in the expansion 
( p6| ) affects only the expression connecting the surface 
impedance with the real surface profile, the scattering 
formulation being rigorous and energy conserving (recall 
that losses are not accounted for) whatever the surface 
impedance is. Nonetheless, inasmuch as we wish to be 
able to quantitatively relate our results with real defect 
sizes, the effect of neglecting the higher order terms in 
Eq. ( p6| ) has to be determined. We have thus verified 
in the most unfavorable cases that including the first or- 
der term in [—id{uj)QY' in Eq. ( p7[ ) barely modifies our 
calculations. 

In order to establish the accuracy and efficiency of the 
numerical calculations based on the formulation above 
we first calculate the function T{q,uj) [cf. Eq. (| 
following the numerical procedure outlined in Sec 
and the Appendix, for two Gaussian defects of half- width 
a/ A = 0.1 and heights h/\ ^ ±0.05 (protuberance and 
indentation of equal height /depth) , where A is the wave- 
length of the SPP. From these results, the SPP reflection 
and transmission coefficients are straightforwardly cal- 
culated [cf. Eqs. (^) and (|24|)], along with the DRC 
[cf. Eqs. ([1^) and (pi|)]. Furthermore, the magnetic and 
electric fields at any point in the vacuum half-space can 
be calculated from Eqs. (|l|) and (|6|) by using Eq. (10). 
In Fig. H we present the results thus obtained for the 
magnetic field intensities at the vacuum-metal interface 
in the vicinity of the Gaussian defects, and for the an- 
gular distribution of the scattered field in the far field. 
From the surface magnetic field in Fig. ||(a), it is evident 
that both surface defects reflect back part of the incom- 
ing SPP, which interferes with the incoming SPP giving 
rise to the oscillatory pattern to the left of the defect 
(negative xi-axis). Near the defect the magnetic field is 
perturbed. The outgoing transmitted SPP is seen to the 
right of the defect. The SPP reflection and transmission 
coefficients are: i?sp = 0.0025 and Tsp ^ 0.9825 for the 
protuberance, and Rsp = 0.0041 and Tsp = 0.9728 for 
the indentation. In Fig. |^(b) a fairly uniform angular dis- 
tribution of the DRC is observed (this will be discussed 
below). The total scattered power calculated from Eq. 
(|^) is S* = 0.0149 for the protuberance, and S = 0.0231 
for the indentation. Energy conservation is thus satisfied 
within a 0.01% error. 

We find that away from the vicinity of the defect, the 
magnetic field is fully described by either the interference 
between incoming and refiected SPP on the left-hand 
side, or by merely the transmitted SPP on the right- 
hand side. This corroborates, as expected, our argument 
in Sec. [ID leading to Eqs. (p^). The reflection and 



A. Energy balance dependence on defect size 

The question now arises naturally as to how efficient 
the surface defect is in coupling the incoming SPP into 
the different outgoing channels (either SPP or photons), 
or conversely, what the appropriate defect parameters 
are that maximize or minimize those channels; this is 
crucial for both an understanding of the scattering pro- 
cess and the design of practical devices. To that end, 
we have studied the dependence of the scattering coef- 
ficients Rsp, Tsp, and S on the defect half- width a for 
both Gaussian protuberances and indentations of differ- 
ent heights h/X = 0.05 and 0.2. The results are shown 
in Fig. 1^. Several general features are evident from these 
results. 

First, SPP reflection is relevant only for very nar- 
row defects, a < A/5, for either protuberances or in- 
dentations. Indeed there iS|-an optimum defect width 



,ximum.a These defects are called 



for which Rsp i: 
plasmon mirrors}^^!^ For increasing defect widths, pro- 
tuberances and indentations begin to behave differently, 
except for their negligible contribution to SPP reflection. 
On the one hand, SPP transmission through protuber- 
ances monotonically diminishes at the expense of radi- 
ation. The conversion is steeper the higher the defect 
is. Indentations, however, exhibit an oscillatory pattern 
with increasing defect width, in such a way that radi- 
ation (SPP transmission) increases (decreases), passes 
through a maximum (minimum) , and then tends asymp- 
totically to (1). The oscillation period, the defect width 
that yields maximum radiation, and the value of this 
maximum, all depend on the surface height. Note that 
both protuberances and indentation may behave as light- 
emitter£3 (high SPP-light conversion efficiency) for an 
appropriate (and distinct) range of defect parameters. 
Below we analyze in detail the behavior of SPP mirrors 
and light-emitters. 



B. Narrow defects: SPP mirrors 

Surface defects playing the role of SPP mirrors hajis 
been studied experimentally in PSTM configurations.Ej 
This phenomenon has been analyzed for different defect 
shapes in Rcf. |3^, where in addition a simple analyti- 
cal prediction is given through a perturbation-theoretic 
argument. In the case of Gaussian-shaped defects, the 
predicted half-width that yields maximum refiection is: 
Gmirr ~ [2^^'^ k-^{u))]~^ . Our numerical results further 
corroborate this prediction, since it is seen in Fig. I that 
a-mirr/^ ~ 0.1 no matter what the defect height is [as 
long as this height does not exceed the range of validity 
of Eq. (psi)]. Nonetheless, the maximum SPP reflection 
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increases with the defect height, and is shghtly larger for 
indentations. 

The electric and magnetic near-field intensities for 
Gaussian defects, placed at the origin, of width a/A = 0.1 
and heights h/X — ±0.05 (protuberance and indentation, 
respectively) are presented in Figs. ^ and ^, respectively,. 
The near-field maps are quite similar in both cases. The 
oscillations to the left of the obstacles clearly reveal the 
interference between the incident and backscattered SPP, 
their period being T « 2TT/{2k^), as expected, and their 
contrast being related to Rsp- A bright region is seen 
to the right of the defects that is due to the strong SPP 
transmission. Poynting vector maps superimposed on the 
electric near-field intensity maps confirm the description 
of the energy fiow given above. 

The corresponding angular distributions of scattered 
light (DRC) have been shown in Fig. ||(b). There are 
no significant qualitative differences between protuber- 
ances and indentations, both yielding a fairly structure- 
less angular dependence; quantitatively, an indentation 
leads to stronger light coupling. The qualitative be- 
havior is somewhat expected: the same perturbation- 
theoretic argument predicting maximum SPP reflection 
for a defect width that maximizes tie scattering poten- 
tial [cf. Eq. (||)] at backscatteringjEj leads to a mostly 
uniform SPP coupling into EM waves in the radiative 
region {\q\ < lu/c). 

C. Wide protuberances: Total light-emitters 

The ability of Gaussian-shaped protuberances to cou- 
ple SPP into light has been pointed out in Ref. Here 
we analyze in detail the conditions for protuberances and 
indentations alike to behave as light-emitters with cou- 
pling efficiencies beyond 90%, larger than that reported 
in Ref. |4[ Figure ^ above illustrates the discussion. 

To begin with, let us focus on protuberances. For 
widths beyond those producing significant SPP reflec- 
tion (SPP mirrors), SPP-light conversion increases mono- 
tonically whereas, as expected from energy conservation, 
SPP transmission decays. This variation is faster for 
higher protuberances. Indeed, the curves in Fig. I (bot- 
tom) indicate that lima^oo S — 1 even for the small pro- 
tuberance. We have found coupling efficiencies beyond 
90% in the case of h/X = 0.2 and a/A > 3.6. In Fig. |, 
the electric and magnetic near-field intensity maps for 
one such defect are shown. The absence of oscillations 
to the left of the defect reveals that SPP reflection is 
small; SPP transmission is small too (though consider- 
ably larger than Rsp), as seen on the right-hand side of 
the defect. A light beam is observed leaving the surface 
from the defect at near-grazing scattering angles. This 
energy flow picture is further corroborated by the an- 
gular distribution of the DRC shown in Fig. 0, with a 
maximum at 9s ~ 74°. Qualitatively, the fact that the 
metal protuberance enters the vacuum half-space seems 



to favor the SPP-photon coupling (playing the role of a 
launching platform) . 

D. Wide indentations: Light-emitters and SPP total 
transmission 

In the case of wide indentations, however, the behav- 
ior of the different outgoing channels differs from that 
for protuberances, and exhibits a richer phenomenology. 
Upon increasing the width of the indentation beyond the 
range of signiflcant SPP reflection (see Fig. ^) , SPP trans- 
mission reaches a minimum value leading to maximum 
radiation, and then slowly grows towards total transmis- 
sion (no radiation) in an oscillatory manner. The defect 
width that yields maximum radiation, its value, and the 
oscillations depend on the defect height. 

To understand such behavior, we plot in Figs. H and ^ 
the electric near-field intensity maps for the higher de- 
fects {h/X = 0.2) of widths a/A = 0.3 and a/A = 4, 
respectively. These widths correspond to the first abso- 
lute and sixth subsidiary, maxima in Fig. ^ (bottom), 
respectively. Both indentations give rise to a negligible 
amount of SPP reflection (no oscillations to the left of the 
defect), as expected from Fig. I (top). SPP tr ansmission 
(to the right of the defect) is very small for a/A = 0.3, 
but a strong light beam at grazing scattering is observed 
{S = 94%). For a/ A = 4, although most of the en- 
ergy goes into Tsp = 86.1%, a small amount of radia- 
tion also at grazing scattering angles is seen [recall that 
even though for this width a local maximum occurs in 
S, its value is very small S — 13.3%, see Fig I (bot- 
tom)]. But what is very illustrative to the discussion on 
the behavior of the outgoing channels is the near-fleld 
within the indentation (strictly speaking, right on top of 
the indentation region, since we are using an IBC on a 
flat surface). Oscillations are found therein, the number 
of minima (one in Fig. ^ and six in Fig. ^ being directly 
related to the position of the corresponding maxima in 
the S vs. a curve [see Fig. ^. Therefore, it can be in- 
ferred that the oscillatory behavior of SPP transmission 
and conversion into light in indentations is governed by a 
cavity-like effect. In fact, the near field map (not shown 
here) in the vicinity of any minimum-radiation indenta- 
tion provides further evidence for this suggestion. 

As a consequence, the range of defect widths for which 
high coupling efficiencies are encountered is far more re- 
strictive for indentations. In fact, in contrast with protu- 
berances, only sufficiently deep indentations {h/X > 0.2) 
are capable of producing radiation efficiencies S > 90%, 
and only for a narrow range of parameters. It seems as 
if the indentation geometry would somehow hinder graz- 
ing light scattering. Therefore, this kind of light- emitters 
might not correlate with any of the reciprocaJr-versions 
(SPP flashlights) seen in PSTM experiments 
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E. Large width limit 

Although the energy conservation criterion is reason- 
ably well satisfied in our calculations, even for defects 
wider than those used in Fig. || (we have reached up to 
a/A = 20), one has to be careful when interpreting the 
results in the limit a/A —>■ oo. It turns out that the 
determination of the behavior of defects in this limit is 
important, since different tendencies have been encoun- 
tered for protuberances and indentations (total radiation 
and transmission, respectively). 

The analysis of the appropriate defect width that yields 
maximum coupling is not simple even if making use of 
the Born approximation, since it requires the evaluation 
of the integral of s{q — k^) for all homogeneous waves 
\q\ < uj/c. And yet such an approximation does not prop- 
erly describe the formally exact numerical calculations. 
Alternatively, we have carried out an analytic calcula- 
tion based on the use of a boundary condition similar to 
the Kirchhoff approximation. The approach relies on the 
expression for the scattering amplitude in terms of an in- 
tegral equation along the surface with the magnetic field 
and its normal derivative inside the integrand (cf. Refs. 
|37| , ^ for the integral equation formulation, and Ref. ^ 
for its version making use of the IBC on a flat surface). 
By assuming that the surface magnetic field is given by 
the incoming SPP, the scattering amplitude reads in the 
large width limit: 



lim R{q,Lu) 

a/A — >C30 



2TrS{q - fc^(w))(l + 7rso)- 



(30) 



Although it does not satisfy energy conservation (not sur- 
prisingly, due to the approximation involved) , the former 
result gives an estimation of the limiting behavior shown 
above (see Fig. p|) : SPP transmission saturates for inden- 
tations (so > 0), whereas protuberances (sq < 0) tend 
to decrease SPP transmission (in agreement with the nu- 
merical calculations, yet the exact limit is not predicted). 



IV. CONCLUSIONS 

We have presented a theoretical formulation that de- 
scribes in a rigorous manner the scattering of a surface 
plasmon polariton propagating along a planar vacuum- 
metal interface by a one-dimensional obstacle modeled 
through an impedance boundary condition. By solv- 
ing the fc-space scattering integral equation upon which 
the formulation is based, the angular spectrum of the 
scattered electromagnetic field in the vacuum half-space 
above the metal surface can be calculated, which in turn 
allows us to obtain the near electric and magnetic fields, 
the amplitudes of the refiected and transmitted SPP, and 
the angular distribution of the intensity of radiated waves 
resulting from the conversion of SPP into volume waves. 
A numerical method to solve the scattering integral equa- 
tion has been put forth. 



We have made use of these calculation methods to 
study the SPP scattering by one-dimensional Gaussian 
defects, either protuberances or indentations. In partic- 
ular, the dependence of the scattering process on the sur- 
face defect parameters has been analyzed. Several con- 
clusions can be drawn from our results with respect to 
the behavior of Gaussian protuberances or indentations. 

SPP reflection is only significant for very narrow sur- 
face defects, with half- widths a < c/uj. Our near field 
results explicitly show that in this case protuberances 
and indentations behave alike, the latter reflecting SPP 
slightly more efficiently. The dependence of the SPP re- 
fiection coefficient on the half-width confirms for different 
defect heights the condition predicted in Ref. ^ of maxi- 
mum SPP reflection, ICjading to the plasmon mirrors seen 
in PSTM experiments. E3 

For wider Gaussian defects, protuberances and inden- 
tations yield a entirely different picture, the only common 
feature being the negligible contribution to SPP reffec- 
tion. Protuberances, on the one hand, increasingly radi- 
ate more light at near grazing scattering angles at the ex- 
pense of SPP transmission. They behave as light- emitters 
with coupling efficiencies approaching 100% with increas- 
ing half- width. The higher the defect is, the larger the 
SPP-light conversion. On the other hand, indentations 
tend to total SPP transmission without radiation with 
increasing half-width. The increase (decrease) of the 
SPP transmission (light coupling) occurs in an oscilla- 
tory manner starting from an absolute minimum (max- 
imum) in transmission (radiation) for small half-widths, 
the period of the oscillations being related to the defect 
impedance in a way reminiscent of a cavity-like effect. In- 
terestingly, we have found that for sufficiently deep inden- 
tations, this maximum radiation value can be extremely 
large (even larger than 90%), so that the Gaussian in- 
dentation thus behave as a light-emitter. 

Our results and discussion provide a thorough picture 
of the different aspects of SPP scattering by surface de- 
fects which, besides being interesting in itself as a scatter- 
ing prrir'ogg appparQ tn hp ug'^fnl in a number of related 
problems.00BEll'mSE[S It is rigorous for one- 
dimensional defects and indeed sheds light on the two- 
dimensional case, and can in turn explain and predict ex- 
perimental results.E^l2j In this regard, it would be inter- 
esting to perform experiments on metal surfaces with de- 
fects of controled profile. With respect to the 2D case, it 
shoukLbe emphasized that the recent work by Shchegrov 
et alEj for circularly symmetric surface defects repro- 
duces the radiation pattern with peculiar lobes in t 
azimuthal angle dependence observed experimentally. 
However, further theoretical work is needed that could 
address more complicated geometries used in the exper- 
iments and/or unexnlaiped processes involving surface 
plasmon-polaritons.liaEaO 
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APPENDIX 

By discretizing q and p in Eq. ( [Tl| ) and replacing the 
infinite limits in the integral by sufficiently large finite 
limits, the following system of linear equations is ob- 
tained for the tn = T{pn,Lu) unknowns: 



(Al) 



with Vm = V {qm\k^ {(^)) ■ The matrix elements Kmn are 
given by 



mn ; 



(A2) 



where 6, 

M, 



mn is the Kronecker delta and 
1 



C(±k^{Lu),L0)V(qrn\±k^{uj)) 

2k^{u}) - Aq/2 



27r 

X ^7ri + log 
1 



2fc«(a;) + Aq/2 



(A3a) 



= -7;-C'(g„i,a;)F(g„i|p„) 

p„ 4- Ag/2 - k'^iuj) 



X l^log 
+ log 



Pn - A(7/2 - k^iiu) 
Pn - Aq/2 + k^iuj 



Pn + Aq/2 + kR{Lj) 

p„ ^ ±fc«(c^). 



(A3b) 



In obtaining Eqs. (A3), the explicit form of Go{q,Lo) 
shown in Eq. (|l^) has been taken into account in cal- 
culating the corresponding integrals over the sampling 
intervals [p„ — Ag/2,p„ -f- Ag/2]. 
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FIG. 1. Illustration of the scattering geometry. 
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FIG. 2. (a) Square modulus of the total surface magnetic 
field and (b) DRC resulting from the scattering of a SPP of 
frequency Tiuj =1.96 eV (A =632.8 nm) by Gaussian defects 
on a silver surface (e = —17.2) of half-width a = O.IA. Solid 
curve: h = 0.05A (protuberance); dashed curve: h = — 0.05A 
(indentation) . 



C(0 



0.08 

0.04 

0.00 
0.8 

" 0.4 

0.0 
0.8 

0.4 

0.0 



10 



10"' 



10 















f 

1 





- - h/X=0.2 

— h/>i=-0.2 
-- h/X=0.05 

- h/X=-0.05 



0.0 



0.4 




1 



a/?i 



FIG. 3. Normalized SPP reflection and transmission co- 
efficients {Rsp and Tsp, respectively), and total radiated 
energy S, as functions of the Gaussian defect half-width: 
Tiu) =1.96 eV (A =632.8 nm) and e = -17.2. Long-dashed 
curve: h = 0.2A; solid curve: h = — 0.2A; dot-dashed curve: 
h = 0.05A; dashed curve: h = — 0.05A. The inset zooms in 
the reflection coefficient in a semi-log scale for narrow defects. 
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FIG. 4. Near field intensity distribution (in a log scale) re- 
sulting from the scattering of a SPP of frequency JvjJ =1.96 eV 
(A =632.8 nm) by a Gaussian protuberance on a silver surface 
[e = -17.2) of half-width a = O.IA and height h = 0.05A. (a) 
Magnetic field intensity and (b) Electric field intensity and 
Poynting vector. 
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FIG. 5. Same as in Fig. W but for a Gaussian indentation 
h = -0.05A. 




-4 4 

X (micronsL 

FIG. 6. Same as in Fig. ^ but for h = 2\ and a — AX (and 
for a larger near field area). 
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FIG. 8. Same as in Fig. |(b) but for h = -2A and a = 0.3A. 
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FIG. 9. Same as in Fig. |(b) but for h = -2A and a = 4A. 
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FIG. 7. DRC as in Fig. g(b) but for: 
h = 2X and a — 4A; dashed curve, h = 
lower solid curve, h = — 2A and a — 4A. 
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